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MAGNETOSHEATH  EFFECTS  ON  CYLINDRICAL 
LANGMUIR  PROBES 


I.  INTRODUCTION 

The  Langmuir  probe  is  one  of  the  most  widely  used  diagnostic  tools  in  laboratory  and 
space  plasmas,  largely  because  of  its  apparent  simplicity  and  the  extensive  documentation 
of  its  operational  characteristics.  However,  its  apparent  operational  simplicity  belies  poten¬ 
tial  pitfalls  which  can  seriously  degrade  its  diagnostic  capability  and  associated  determina¬ 
tions  of  plasma  density  and  temperature1'7.  Some  appreciation  for  potential  problems  can 
be  achieved  by  recognizing  that  the  current-voltage  (I-V)  characteristic  of  the  probe  is 
a  function  of  the  size,  geometry,  and  surface  condition  of  the  probe,  as  well  as  the  probe 
voltage,  the  ambient  magnetic  field,  and  the  plasma  properties  of  the  charged-particle 
number  densities,  particle  distribution  functions,  and  collision  frequencies. 

The  determination  of  density  and  temperature  from  the  probe’s  I-V  characteristic 
hinges  on  the  experimenter’s  ability  to  accurately  account  for  the  varying  parametric 
influences  and  to  overcome  such  experimental  difficulties  as  plasma  depletion,  contamina¬ 
tion,  and  area  effects.  If  one  assumes  that  all  experimental  sources  of  error  have  been 
eliminated,  the  accuracy  of  the  technique  is  limited  by  the  theoretical  description  of  the 
probe’s  response  under  the  prevailing  plasma  conditions.  Under  these  circumstances,  the 
most  difficulty  is  encountered  when  the  operation  of  the  probe  is  in  any  one  of  the  various 
transition  regions,  where  we  use  the  term  “transition  region”  to  describe  any  domain 
between  the  mathematically  convenient  limits  of  collision-free  and  collision-dominated, 
thin-sheath  and  thick-sheath,  and  strong-field  and  weak-field.  These  regions  are  particularly 
difficult  to  describe  because  one  must  account  for  detailed  charged -particle  trajectories  that 
have  no  convenient  mathematical  form  as  they  traverse  the  region  between  the  undisturbed 
plasma  volume  and  the  probe  surface.  In  this  work  we  look  at  the  transition-region  of 
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magnetic  field  effects  and  the  associated  response  of  cylindrical  Langmuir  probe  electron- 
saturation  currents. 

Probe  response  in  magnetoplasmas  can  be  grouped  into  three  broad  categories  defined 
by  the  relative  magnitudes  of  the  probe  radius  Rp,  the  sheath  thickness  (Rs-Rp),  and  the 
Larmor  radii  for  electrons  (R^)  and  ions  (R‘L).  We  define  these  categories  as: 


V(VV^ReL‘ 

(weak-field) 

(la) 

Rp.<R.-Rp»R£i 

(strong-field). 

(lb) 

(transition-field). 

(lc) 

Each  of  these  categories  has  its  own  morphological  subdivision  established  by  the  independ¬ 
ent  ratios  R^’  ‘/Rp  and  R^’  ‘/(Rg-Rp).  The  first  ratio,  R^’  '/Rp,  involves  geometrical  effects 
which  result  in  magnetic  field  shadowing8  and  the  associated  perturbation  of  a  fully- 
Maxwellian  plasma  distribution  at  the  sheath  edge.  A  number  of  authors2,  912  have  used 
this  ratio  to  describe  magnetic  field  effects  on  cylindrical  probes  and  have  shown  reduced 
saturation  currents  when  the  ratio  was  small.  Miller11  and  Laframboise  and  Rubinstein12, 
however,  infer  that  magnetic  field  effects  may  occur  even  when  R£  1  >  Rp  if  R^’  '  <  (Rg-  Rp ). 
It  is  the  latter  inequality  which  is  of  primary  concern  in  the  present  investigation. 

Arguments  advanced  in  subsequent  sections  point  out  that  the  interrelationships 
between  Rp,  R^>  ',  and  (Rg-Rp)  are  all  important  in  establishing  the  probe’s  response  in 
the  presence  of  magnetic  fields.  Our  results  demonstrate  the  important  role  of  sheath  size 
in  the  case  of  moderate-to-weak  magnetic  fields  and  point  to  possible  errors  that  can  arise 
in  attempts  to  monitor  relative  electron  (ion)  densities  by  changes  in  electron-  (ion-)  satura¬ 
tion  currents.  We  specifically  show  that,  in  general,  one  cannot  assume  le  «  Ne  when  probe 
operation  is  in  the  transition  region  of  magnetic  field  effects. 
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II.  SHEATH  CONSIDERATIONS 

A.  Sheath  size.  In  describing  magnetic  field  effects  on  probes  it  is  necessary  to  establish 
an  understanding  of  the  plasma  sheath  which  envelopes  the  probe  and  controls  the  charged- 
particle  fluxes  as  they  move  from  the  undisturbed  plasma  volume  to  the  collecting  surface 
of  the  probe.  Typically,  the  electron  Debye  length  is  used  as  a  convenient  approximation 
for  the  sheath  size,  but,  this  approach  has  limited  accuracy  and  does  not  include  sheath 
dependence  on  probe  potentials.  We  will  therefore  extract  from  the  numerical  calculations 
of  Laframboise13  an  analytical  relationship  for  the  plasma  sheath  size  which  will  be  applied 
to  our  model  of  magnetic-field  effects  on  cylindrical  probes. 


The  Laframboise  calculations  are  based  on  the  Boltzmann-Vlasov  equation  for 
Maxwellian  electrons  and  ions,  coupled  with  Poisson’s  equation  and  the  boundary  condition 


d^'i 

dR 


1 

R 


(2a) 


for  the  solution  in  the  case  of  cylindrical  probes.  (Additional  discussions  of  Laframboise’s 
calculations  can  be  found  in  the  works  of  de  Leeuw14 ,  Sonin15,  and  Szuszczewicz5.) 

For  electron-attracting  probe  potentials  in  the  limit  of  zero-temperature  ions  the  boundary 
condition  is  relaxed  and  replaced  by 

V  V 

dsfi 
dR 


=  $ 


=  0. 


(2b) 


In  this  case  Rh  is  a  sharp  boundary  defined  as  the  sheath  edge  radius  where  the  repelled  ion 
spatial  distribution  is  discontinuous  [i.e.,  (Nj )R  <-  R  =  Oand(Nj)R>R  =Ne].Inthe 

b  b 

computational  scheme  the  sheath  radius  was  not  known  in  advance,  but  its  position  was 
adjusted  until  the  total  space  charge  within  the  sheath  exactly  canceled  the  charge  on  the 
probe.  The  resulting  boundary  value  problem  then  required  solutions  for  the  potential  and 
charge  density  distributions  as  well  as  the  position  of  the  sheath  edge  Rfa(  =  Rg).  The 
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analytic  expression  which  we  have  fit  to  the  resulting  Laframboise  calculations  (Fig.  38, 
Ref.  13)  for  the  dimensionless  sheath  thickness  S  is  given  by 


R„  -  Rn 
s  p 


2.50-  1.54  exp(-0.32R  /XD) 


\  1/2 


kT 


(3) 


This  function  is  plotted  in  Fig.  1.  Comparison  of  Eq.  (3)  with  the  calculations  of 
Laframboise13  showed  agreement  to  be  within  a  few  percent  for  RpMD  ^  10.  The 
difference  increases  to  26%  at  e<pp/kTe  =  25  in  the  limit  of  Rp/AD  °°;  but  this  deviation  is 
inconsequential. 


Figure  1  also  includes  calculations  for  sheath  radius  based  on  the  analysis  of  Bettinger 

and  Walker16  for  cylindrical  electron-attracting  probes.  Their  analysis  does  not  provide  the 

potential  distribution  within  the  sheath.  Instead,  they  consider  the  probe  surface  and  sheath 

edge  as  concentric  electrodes  of  a  cylindrical  capacitor  and  assume  the  charge  density  within 

the  sheath  to  be  inversely  proportional  to  the  local  velocity  of  the  incoming  particles.  The 

probe  surface  charge  is  then  equated  to  the  total  charge  within  the  sheath  with  the  resulting 

expression  demonstrating  the  dependence  of  sheath  size  on  probe  radius,  Debye  length,  and 

applied  probe  potential.  Although  their  results  do  not  bear  the  rigor  of  the  Laframboise 

calculations,  they  do  provide  an  alternate  description  of  the  sheath  edge  radius  that  includes 

dependence  on  probe  potential  and  radius  rather  than  a  simple  identification  with  Debye 

length.  (It  should  be  noted  that  Eq.  (14)  in  the  work  of  Bettinger  and  Walker16  is  in  error 

by  a  constant,  (s2  +  2ps)ln(l  +  s/p)  =  tt  i//s(1  +  |  4*s)112  should  be  (s2  +  2ps)ln 
2  1/2 

(1  +  s/p)  =  2^(1  +  ■=■  i g)  .  The  corrected  form  has  been  utilized  in  the  current  work.) 

Of  the  two  approaches,  the  Laframboise  results  represent  a  more  realistic  physical 
model  with  a  solution  scheme  that  includes  the  dependence  of  sheath  potential  distribution 
and  the  associated  charge  profile  on  the  geometry  of  the  probe.  It  is  the  Laframboise 
approach  with  the  analytical  description  provided  by  Eq.  (3)  which  is  utilized  in  this  work. 


B.  Magnetic  field  effects.  No  attempt  will  be  made  to  establish  an  anisotropic  description 
of  sheath  geometry,  an  effort  which  goes  beyond  the  scope  of  this  work.  In  this  section  we 


aY/(da  -  sy) 
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is  the  electron  plasma  frequency,  and 


is  the  electron  gyrofrequency. 


cu®  =  eB/mec 


The  criteria  established  in  (1),  particularly  as  they  apply  to  the  effect  of  R?  /( H  -  R  ) 

Li  S  p 

on  electron  saturation  currents,  can  now  be  expressed  as 


p  c 


F^/kT,)1'2 


where 


>  1 

( weak  -magnetosheath ) 

(7a) 

« 1 

( strong-magnetosheath ) 

(7b) 

» 1 

(transition -magnetosheath ) 

(7c) 

1.54  exp  (-0.32  Rp/XD). 

(8) 

These  inequalities  show  that  a  weak-,  transition-,  or  strong-field  classification  depends 
not  only  on  the  magnitude  of  the  magnetic  field  but  also  on  plasma  parameters  of  density  and 
temperature,  as  well  as  the  probe  size  and  applied  potential  <p  .  Thus  a  0.25  G  field  could  have 
a  similar  effect  on  an  ionospheric  plasma  sheath  (typically,  N™14*  =  106  cm'3  and  Te  =  2000°K 
at  F-  region  altitudes)  as  a  30  kG  field  in  a  confined  hot,  dense  plasma  [Np  =  5(1015)cm~3, 

Tp  =  1.16(107  )°K  (=  1  keV)].  Examining  these  conditions  in  terms  of  (7)  we  find 

upK 


F^p/kT/'2 


in  the  ionospheric  case,  and 


U)e  / GJe 

- — - - =  0.95  (10) 

F^p/kT^1'2 

for  the  hot,  dense  plasma.  These  results  assume  e^p/kTe  =  10  as  a  nominal  operational  value 
for  the  collection  of  electron  saturation  currents  by  a  cylindrical  probe  of  radius  Rp  = 

3.8(10  2)  cm.  Since  F,qs.  (9)  and  (10)  are  of  comparable  magnitude,  their  results  show  that 
dramatically  different  probe-plasma  systems  can  have  similar  classifications  with  regard  to 
magnetosheath  effects. 
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C.  Angular  dependence.  The  anisotropic  nature  of  charged-particle  motion  makes  it 
necessary  to  consider  the  field  direction  B  relative  to  the  probe  and  sheath  axes  L.  The  most 
complete  work  done  to  date  in  this  area  is  that  of  Laframboise  and  Rubinstein12  who  have 
conducted  a  theoretical  analysis  of  a  cylindrical  probe  in  a  collisionless  plasma,  with  the 
probe  operating  under  thick  sheath  conditions  at  an  arbitrary  angle  0  £=  cos'1  (  L-B/LB^j 
with  respect  to  a  uniform  magnetic  field.  For  a  probe  at  plasma  potential,  their  analysis  is 
exact;  but  in  regions  of  electron  saturation  currents  their  theory  provides  only  an  upper 
bound  and  an  adiabatic  limit.  These  limits  are  approached,  respectively,  at  larger  and 
smaller  values  of  R?  / L  .,  where  L  .  is  defined  by  Laframboise  and  Rubinstein  as  the  gradient 
scale  length  of  sheath  potential  #/\'Vp\.  The  adiabatic  limit  corresponds  to  R^/L^  O  and 
can  be  represented  approximately  by  R^/(Rs~Rp)  -*•  O. 

The  Laframboise  and  Rubinstein  adiabatic  calculations  show  that  Ic(0a )/Ie(^b )  <- 1  f°r 
9  <  0.  .  A  sample  of  their  electron-current  calculations  at  plasma  potential  (ey? /kT  =  0) 
and  in  the  saturation  regime  (e^p/kTe  =  5)  is  presented  in  Fig.  2.  All  results  are  normalized 
to  fully  Maxwellian  current  levels  collected  at  plasma  potential  under  zero  magnetic  field 
conditions.  In  the  next  section  we  present  our  experimental  results  for  comparison  with 
theory  and  discuss  the  role  of  R^/(Rs-  Rp)  in  approaching  the  adiabatic  limit. 

III.  EXPERIMENTAL  RESULTS 

A.  Experimental  configuration.  The  experimental  results  were  collected  by  a  pulsed- 
plasma-probe  (P3 )  flown  on  a  scientific  rocket  payload  (NASA  18.170)  into  the  equatorial 
ionosphere  above  Lima,  Peru.  The  P3  technique  utilizes  a  special  pulse  procedure  for 
generating  Langmuir-probe  current-voltage  characteristics.  The  procedure  reduces  possible 
perturbations  of  the  measurement  technique  by  surface  contamination6, 17  and  makes 
possible  a  determination  of  electron  density  and  temperature  under  fluctuating  plasma 
conditions18'20. 

Current-voltage  characteristics  were  generated  by  voltage  pulses  (=  100  psec)  which 
followed  a  sawtooth  envelope  from  -1.6  V  to  +  3.0  V.  During  the  interpulse  period 
(a>  1.6  msec)  the  probe  potential  was  held  at  a  fixed  baseline  value  V0  in  order  to  stabilize 
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»  1  p  C'K  p  » 

5T  ^7)  =  ,a51 
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iTfT  =  5 


/  RP  e<t)P  V 

_ ReL  '  kTe 

(0,5)  (0.5,0) _ _ 

(ADIABATIC  LIMIT) _ — - 


_ I  I _ I  I  I 

10°  20°  30°  40°  50°  60°  70°  8 

6  [  =  COS'1  (L  B)] 

Fig.  2  —  Effect  of  probe  orientation  on  electron  current  collection  in  a  magnetic  field 
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probe-surface  conditions.  Changes  n  probe  current  which  might  o(  cur  independent  of 
sweep  voltage  were  tracked  by  a  current  measurement  JB  during  the  time  of  the  blanked 
fifth  pulse6,  1 '  With  VB  =  +  2.QV,  IB  collection  took  place  in  the  electron  saturation 
portion  of  the  I-V  characteristic.  (In  more  recent  applications21 ,  the  fifth  pulse  has  been 
restored  and  IB  -measurement.1  nave  been  performed  during  every  interpulse  period  at 
V  =  VB.) 

A  redundant  pair  of  probes  was  flown  on  the  Nike-Tomahawk  payload.  Each  probe 
was  constructed  of  a  0.076  cm  diani  tungsten  rod  20.3  cm  in  length.  In  its  deployed 
configuration  the  probe  was  perpendicular  to  the  spin  axis  of  the  payload  and  separated 
from  its  39.4  cm  boom  by  a  coaxial  guard  electrode  that  was  1.3  cm  long  and  0.32  cm  in 
diameter.  The  guard  was  simultaneously  driven  with  the  same  voltage  as  the  probes. 

The  payload  was  launched  from  the  Chilca  Range  in  a  trajectory  that  was  4°  north  of 
due  west.  With  the  range  located  at  the  magnetic  equator  and  the  geomagnetic  field 
inclined  at  4°  east  of  north,  the  payload  trajectory  was  perpendicular  to  the  magnetic  field. 
To  maintain  aspect  stability  the  payload  was  spun  at  a  rate  of  4  rps.  The  spin  axis  of  the 
payload  was  inclined  at  10°  from  local  vertical  with  an  estimated  coning  angle  less  than  5°. 
As  a  result,  each  probe  was  rotated  from  0°(*q  )  to  90J(i0°)  with  respect  to  the  magnetic 
field  twice  during  each  spin  period. 

B.  Results.  To  measure  the  dependence  of  IR(=  I®1 )  modulation  on  ionospheric  plasma 
parameters  we  have  plotted  in  Fig  1  values  of  IB  at  0  =  0°  (*q°)  and  0  =  90°  (±0°).  Using 
the  IB(0  =  90°)  profile  as  the  more  accurate  measure  of  relative  density11  and  establishing 
the  conversion  Ne(cm  1 )  =  1.25  (1011  )IR  (A),  Fig.  3  demonstrates  the  importance  of 
plasma  density  tthrough  its  control  of  sheath  size)  in  determining  the  effect  of  magnetic 
fields  on  electron  current  collection  by  cylindrical  Langmuir  probes.  (The  Np/IB 
proportionality  was  determined  near  apogee  by  conventional  analysis2  of  the  electron- 
saturation  portion,  Ijf1,  of  the  current-voltage  characteristic,  i.e.,  dfl^^/dip  -*  Np.  The 
simultaneous  measurement  of  N(  and  IB,  made  possible  with  the  P3  technique,  yielded 
the  constant.  Sources  of  error  identified  with  possible  plasma  depletion7,  surface 


9 


SZUSZCZEWICZ  AND  TAKACS 


I 


contamination17'18,  reference  electrode  area5,  and  convective  effects1  were  inconsequen¬ 
tial.)  In  the  ionospheric  E-region  trough  (125  -  170  km)  where  the  plasma  density  is 
lowest  Qfc  6.0(103 )  cm-3]  the  percent  modulation,  M  =  j[IB(90°)-  IB(0° )]  /IB(90°)|  X 
100%  =  75%,  is  much  greater  than  in  the  F  region  (Z  >  170  km)  where  the  modulation 
is  only  10-15%.  The  difference  is  attributed  to  sheath  size  variation  since  over  the  altitude 
range  in  this  investigation  the  earth’s  magnetic  field  and  associated  Larmor  radii  are 
approximately  constant.  (RL  values  are  constant  only  if  temperatures  are  constant,  a 
situation  which  is  not  generally  true  over  this  altitude  range  where  we  can  expect  up  to  a 
factor  of  two  difference.  But  we  can  neglect  the  temperature  effect  (2X)  compared  to  the 
density  effect  (X100).) 
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The  results  in  Fig.  3  identify  a  potential  problem  area  for  plasma  experimenters  who 
utilize  fixed  bias  cylindrical  probe  measurements  of  electron -saturation  currents  to 
determine  changes  in  electron  density.  Even  when  the  probe  is  held  at  a  fixed  angle  with 
respect  to  the  magnetic  field,  the  spatial  or  temporal  profile  of  plasma  density  can  be 
distorted  by  changing  sheath  sizes  that  accompany  varying  plasma  densities.  Distorted  data 
can  result  in  misleading  interpretations  of  active  physical  pnnciples.  In  Fig.  3  the  Ifi(0  =  0°) 
curve  could  lead  to  erroneous  conclusions  concerning  nighttime  E-region  depletion 
mechanisms  (130-170  km)  or  applicability  of  the  electron  density  gradient  at  the  bottom- 
side  of  the  F  layer  (170-240  km)  to  the  Rayleigh-Taylor  instability  and  the  triggering  of 
ionospheric  plasma  irregularities. 

That  reasonable  levels  of  measurement  integrity  can  still  be  maintained  in  weak 
magnetoplasmas  has  been  indicated  by  Miller11,  and  Laframboise  and  Rubinstein12.  Their 
works  indicate  that  accuracy  is  improved  in  cylindrical  probe  measurements  when  the  angle 
between  the  probe  and  the  magnetic  field  is  large  (0  60°).  This  is  particularly  true  in  the 

case  e<£p/kTe  =  0  and  Rp/R^  <  0.5  as  shown  by  the  solid  curves  in  Fig.  2.  (The  cylindrical 
orbit-limited  boundary  is  included  in  the  figure  because  it  is  a  realistic  upper-bound  that 
should  not  be  exceeded.  Laframboise  and  Rubinstein  make  this  point  in  their  discussion.) 

Three  cases  have  been  selected  from  our  ionospheric  data  to  detail  the  behavior  of 
1B(0)  as  a  function  of  plasma  density.  The  results  have  been  normalized  to  1FJ(90°)  and 
plotted  in  Fig.  4  as  curves  A,  B,  and  C.  We  note  that  the  modulation  increases  with 
decreasing  Np ,  a  parametric  dependence  not  shown  in  current  theories  involving  thick 
sheath  conditions  (RpAD  ^  1).  Specifically,  we  find  that  (ojp/u>*)/F  ^e^p/kTp^  1/2 
£  =  R^/(Rs~Rp)J  equals  2.5,  0.47,  and  0.3  for  A,  B,  and  C,  respectively.  In  terms  of  the 
inequalities  in  (7)  these  cases  qualify  as  transition-magnetosheath.  We  observe  that  the 
R^/IR^-Rp)  =  2.5  case  has  the  smallest  modulation  since  it  approaches  the  condition  of 
weak -magnetosheath.  The  data  show  that  the  modulation  would  not  be  zero  as  a  result  of 
Rp/R^  -»  0  alone.  The  modulation  can  be  zero  only  if  Rp/R^  and  (Rs-Rp)/R^  both  go 
to  zero,  a  combined  condition  represented  by  R^/(Rs-Rp)  >  1  (Eq.  (7a)]  in  the  thick 
sheath  limit.  The  data  demonstrate  the  important  coupling  of  B,  Np,  Tp,  and  ifip  in 
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determining  the  degree  to  which  magnetic  fields  («  rturli  electron  current  collection.  One 
cannot  give  sole  consideration  to  H  /Xjj  or  Rp ' K y  ,  hut  rather  their  important  interrelation¬ 
ships  as  described  in  Eq.  (7 ). 

C.  Comparison  with  theory  Figure  4  includes  the  upp.  r  bound  (l1)  and  adiabatic  limit  (Li 
predictions  of  Laframboise  and  Rubinstein  wih  parametrc  dependence*  on  R  /RJ The 
adiabatic  limit  results  are  based  on  two  models  For  non  zero  values  of  R  /R?  the 
calculated  current  results  from  a  detailed  study  of  the  intersections  of  helical  orbits  with  an 
infinitely  long  probe.  Current  is  collected  only  from  positive-total-energy  orbits  that 
originate  at  infinity  rather  than  elsewhere  on  the  probe.  The  particles  are  assumed  to  have  a 
“one-dimensional”  velocity  distribution  in  which  kinetic  energy  gain  from  the  applied  field 
goes  entirely  into  increased  speed  along  B.  In  this  model  all  possible  collection  orbits  are 
empty  and  the  associated  probe  current  is  zero  when  the  probe  axis  is  parallel  to  B.  The 

calculations  of  Laframboise  and  Rubinstein  for  R  /R?  =  0.1,  1.0,  and  10  are  shown  in 

P  L 

Fig.  4  with  extrapolations  to  0  =  0° . 

At  Rp/R‘L  =  0  a  different  approach  is  employed  by  Laframboise  and  Rubinstein.  As 
already  noted,  the  current  is  calculated  in  the  limit  in  which  the  velocity  distribution  of 
incoming  particles  at  the  probe  is  “one-dimensional”;  but,  instead  of  detailed  orbital 
considerations,  their  scheme  is  reduced  to  calculating  the  current  collection  of  particles 
with  orbits  of  negligible  curvature  that  flow  along  magnetic  flux  tubes  intersecting  the 
probe.  The  currents  are  allowed  to  have  components  from  particle  velocities  both 
perpendicular  and  parallel  to  B.  This  is  an  important  difference  since  it  allows  current 
collection  even  with  the  probe  aligned  with  B.  Compared  with  the  calculations  at 
Rp/R*  =£  0,  the  approach  at  Rp/R^,  0  =  0  effectively  eliminates  the  emptying  of  collec¬ 
tion  orbits  by  previous  intersections  with  an  infinitely  long  probe.  The  result,  as  shown  in 
Fig.  4,  is  a  finite  nonzero  value  for  current  when  the  probe  is  parallel  to  B.  In  a  real  situa¬ 
tion  with  a  probe  of  finite  length,  collection  orbits  at  6  =  0°  would  be  partially  filled  by 
particles  entering  the  sheath  at  the  probe  ends.  One  would  therefore  expect  current  collec¬ 
tion  at  small  values  of  0  and  Rp/R^  to  be  bracketed  between  the  predictions  of  the  two 
computational  schemes. 
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Fig.  4  —  Experimental  (curves  A,  B,  and  C)  and  theoretical  (curves  U  and  L  with  parametric  dependence  on 

R  / R*  )  results  showing  the  effect  of  magnetic  aspect  on  cylindrical  electron-collecting  probe  current  The 
p  L 

L -curves  at  0  ^  6°  are  taken  from  Ref.  12  with  dotted  extrapolations  to  0  -  0°  based  on  Eq  (6)  in  the  same 
reference 
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The  experimental  results  in  Fig.  4  (Curves  A,  B,  and  C)  were  all  collected  under  condi¬ 
tions  for  which  R  /R?  *  0.014  and  are  therefore  appropriate  for  comparison  with  the 
Laframboise  and  Rubinstein  curves  U,  and  L(Rp/R^,  6)  =  L(0,  0),  L(0.1,  0).  The  data 
fall  between  the  theoretical  limits  (with  the  qualification  on  small  6  discussed  here)  and 
provide  evidence  of  a  continuous  variation  from  conditions  of  a  weak-magnetosheath 
(the  upper  bound.  II)  to  those  of  a  strong-magnetosheath  (the  adiabatic  limit,  L).  The  trend 
in  the  data  also  supports  the  premise  advanced  in  the  previous  section  that  Rl/(Rs"  Rp) 
is  an  approximate  measure  of  R? /fa/l^^l)  where  <^/  i  Vf  |  is  the  gradient  scale  size  of  sheath 
potential  The  analytic  formulation  in  Eq.  (7)  is  therefore  a  useful  gauge  for  determining  the 
degree  to  which  adiabatic  limit  calculations  apply  to  a  given  experimental  condition. 

The  trend  in  the  data  and  the  minimum  observed  value  of  R[  /(Rg-  Rp)  indicate  that 
the  experimental  conditions  did  not  meet  the  guidelines  for  having  achieved  the  adiabatic 
limit  R|  (Rg-R(i)<  1  We  would  infer  that  with  increasingly  smaller  values  of  R^/(Rg- Rp) 
experimental  curves  operating  under  thick  sheath  conditions  would  approach  L(Rp/R^,  6)  = 
L(0  1.  for  0  £  5  At  0  -  0°.  a  finite  nonzero  current  would  result  such  that  0  <  IB(0°)/ 
IB(90r  )<  0.28 

Finally,  we  note  that  at  no  time  in  our  experiment  did  we  observe  a  “negative 
resistance"  current  voltage  characteristic,  a  feature  predicted  by  the  detailed  orbital  scheme 
of  Laframboise  and  Rubinstein  1  his  supports  their  position  and  our  observation,  namely 
that  the  requirements  for  achieving  the  adiabatic  limit  [our  inequality  (7h)]  are  very  severe. 

IV.  COMMENTS  AND  CONCLUSIONS 

Arguments  have  shown  that  the  extent  to  which  a  magnetic  field  modifies  the  response 
of  a  cylindrical  Langmuir  probe  cannot  be  measured  solely  by  the  magnitude  of  the  field 
itself.  In  the  past,  primary  concern  has  been  given  to  the  ratio  of  probe  radius  to  the  Larmor 
radius  in  efforts  to  define  weak  ,  strong  ,  and  transition-field  response.  We  have  shown  that 
an  equally  important  consideration  is  the  ratio  of  Larmor  radius  to  sheath  thickness,  and 
have  established  appropriate  inequalities  (7a),  (7b),  and  (7c)  involving  not  only  the 
magnetic  field  strength,  but  also  the  plasma  parameters  of  density  and  temperature,  and  the 
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applied  probe  potential.  Our  results  suggest  that  the  regions  of  operation  are  more  appro¬ 
priately  defined  as  weak-,  strong-,  and  transition-magnetosheath  in  order  to  more  clearly 
stipulate  the  interacting  phenomena.  We  have  shown  that  a  0.25  G  field  can  have  the  same 
perturbing  influence  on  a  cylindrical  probe  applied  to  ionospheric  plasma  densities  as  a 
30  kG  field  would  have  in  a  hot,  dense  laboratory  plasma. 

Our  experimental  results  are  in  agreement  with  the  criteria  [Eq.  (7a-c)]  which  define 
the  various  regions  of  response.  In  particular,  the  results  reveal  potential  problems  in  the 
use  of  fixed-bias  probes  as  indicators  of  changes  in  relative  electron  densities  in  a 
magnetoplasma.  Because  sheaths  expand  and  contract  as  plasma  densities  fall  and  rise,  the 
ratio  of  gyroradius-to-sheath  thickness  varies  and  modifies  the  probe  response  in  a  non¬ 
linear  way.  One  must  therefore  guarantee  weak-magnetoplasma  conditions  (7a)  before 
assuming  Ie  =  KNg. 

Our  results  have  also  demonstrated  the  importance  of  sheath  thickness  on  the  response 
of  cylindrical  probes  as  a  function  of  probe  orientation  relative  to  the  ambient  magnetic 
field.  We  have  shown  that  the  percent  modulation,  i.e.,  M  =  100  X  £le(90°)  -  Ie(0°)^  / 
Ie(90°).  increases  with  sheath  thickness.  This  result  is  in  agreement  with  the  upper  bound 
and  adiabatic  limit  calculations  of  Laframboise  and  Rubinstein  and  provides  useful  guide¬ 
lines  for  determining  the  degree  to  which  the  opposing  limits  are  met.  In  particular,  a 
comparison  between  theory  and  experiment  shows  that  our  weak-  and  strong- 
magnetosheath  conditions  are  in  effect  identical  to  their  upper  bound  and  adiabatic  limits 
respectively.  Using  the  criteria  expressed  by  the  inequalities  in  (7),  our  experimental  condi¬ 
tions  did  not  satisfy  the  guidelines  for  having  achieved  the  adiabatic  limit,  a  result  consistent 
with  the  non-existence  of  a  “negative  resistance”  characteristic  (dl“l/d^  <  O)  in  our 
data.  This  supports  the  Laframboise  and  Rubinstein  position  that  the  requirements  for 
achieving  the  adiabatic  limit  [our  inequality  (7b)]  are  very  severe. 

Finally,  our  results  indicate  that  Ig(0  *  90°)  is  a  better  measurement  of  electron 
density  variations  than  Ig(0  <4  60°).  However,  without  an  independent  measurement  of 
Ne,  it  is  impossible  to  establish  the  absolute  accuracy.  We  suggest  that  this  area  be 
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approached  with  caution,  particularly  any  use  of  standard  analysis  techniques  in  the 
electron  saturation  region  for  making  absolute  measurements  of  ambient  electron  density. 
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